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' Abstract 

■ In this paper we examine a new class of five dimensional (5D) 

£SJ . exact solutions in extra dimension gravity possessing Lie algebroid 

symmetry. The constructions provide a motivation for the theory 
of Clifford nonholonomic algebroids elaborated in Ref. pQ. Such 
Einstein-Dirac spacetimes are parametrized by generic off-diagonal 
metrics and nonholonomic frames (vielbeins) with associated nonlin- 
Q H ' ear connection structure. They describe self-consistent propagations 

of (3D) Dirac wave packets in 5D nonholonomically deformed Taub 
NUT spacetimes and have two physically distinct properties: Fist, the 
metrics are with polarizations of constants which may serve as indi- 
^ ■ rect signals for the presence of higher dimensions and/or nontrivial 

torsions and nonholonomic gravitational configurations. Second, such 
Einstein-Dirac solutions are characterized by new type of symmetries 
defined as generalizations of the Lie algebra structure constants to 
nonholonomic Lie algebroid and/or Clifford algebroid structure func- 
tions. 
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1 Introduction 

Recently we studied a new class of solutions in string gravity (with nontrivial 
limits to the Einstein's gravity) possessing Lie algebroid symmetry [2] and 
describing black holes in solitonic backgrounds, see Ref. [3] for details on 
algebroid theory and related bibliography. These solutions were generated 
using the method of anholonomic frames with associated nonlinear connec- 
tion (N-connection) structure |E]. The technique can be extended for spinor 
variables and Einstein-Dirac spaces and motivates the concept of Clifford 
algebrois [T]. In this paper we apply such spinor and Lie algebroid methods 
in five dimensional (5D) gravity in order to construct metrics defining non- 
holonomic deformations of the Taub NUT metric to certain Einstein-Dirac 
algebroids defined by exact solutions of the Einstein-Dirac equations with 
Lie algebroid symmetry. 1 

1 Various classes of 4D metrics induced from 5D Kaluza-Klein theory and brane/string 
gravity are involved in many modern studies of physics related to gravitational instantons 
and monopoles and to the geodesic geometry of higher dimension spacetimes [S] ; the con- 
structions where generalized to generic off-diagonal and locally anisotropic gravitational 
configurations in Refs. [HJ. 
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Exact solutions pay a special role in classical and quantum gravity pro- 
viding a testing ground for fundamental concepts and general methods which 
can be studied by approximation methods. Moreover, any approximations 
are usually based on some exact solutions. An exact solution of the Ein- 
stein equations is characterized by corresponding symmetries and boundary 
(asymptotic) conditions. For instance, in modern astrophysics the asymp- 
totically Minkowski/de Sitter metrics, with spherical or cylindrical symme- 
try, are of special interests. In cosmology, additionally to the spherically 
symmetric solutions, one considers anisotropic models defining spacetimes 
with Lie group symmetry. Nevertheless, following various fundamental pur- 
poses in quantum gravity and string theory, it is important to investigate 
more general spacetime models described by generic off-diagonal metrics 
and nonholonomic structures, deformed (non) commutative symmetries and 
non-perturbative gravitational vacuum configurations (solitons, instantons, 
monopoles, spinor and pp- waves,...). In this line, the approach tp construct- 
ing and investigating classes of solutions with Lie algebroid symmetry dis- 
tinguishes a new direction in mathematical gravity and string theory. Such 
spacetimes preserve a number of important features of manifolds with group 
symmetry but posses a more reach geometric structure combining the prop- 
erties of bundle spaces and various type of nonlinear symmetries, singular 
maps and nonholonomic constraints. 

There are certain applications of the algebroid theory in geometric me- 
chanics [7j and, recently, there were elaborated some algebroid approaches 
in the theory of gauge fields, gravity and strings and noncommutative ge- 
ometry [HI El HI 12]- Perhaps, one can reflect on groupoid and algebroid 
program of geometrizations of physics considered as a modern versions of 
the Felix Klein's "Erlanger Propgram" (1887) when instead of groups and 
algebras one considers, respectively, groupoids and algebroids. In a more 
general context, containing spinors and Clifford algebras, one has to extend 
the constructions to the so called C-space, Clifford space, i.e. the space of 
Clifford numbers, or Clifford aggregates (see details and references in [TU]). 
but provided with additional geometric structures, for instance, of Clifford 
nonholonomic algebroid [2], in order to describe certain classes of nonlinear 
gravitational interactions with generalized symmetries. 

Motivated by the mentioned developments and prospects, in this paper 
we study some explicit examples of exact solutions in gravity when the Lie 
algebroid and Clifford algebroid structures are modelled as Einstein-Dirac 
configurations defined by nonholonomic deformations of Taub NUT metrics. 
The Hawking's j3] suggestion that the Euclidean Taub-NUT metric might 
give rise to the gravitational analogue of the Yang-Mills instanton holds 
true for generalization to "algebroid instantons". Neverthaless, in this case 
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the solutions have some anisotropically polarized constants being of higher 
dimension gravitational vacuum polarization origin. These nonholnomically 
deformed Taub-NUT metrics also satisfy the vacuum Einstein's equations 
with zero cosmological constant when the spherical symmetry is deformed, 
for instance, into an ellipsoidal/ toroidal configuration or transformed into 
locally anisotropic wormhole metrics, see details in Refs. jHj- Such algebroid 
Taub-NUT metrics can be used for generation of deformations of the space 
part of the line element defining Lie algebroid modifications of the Kaluza- 
Klein monopole solutions proposed by Gross and Perry and Sorkin [TT] . 

The Schrodinger quantum modes and the Dirac equation in the Euclidean 
Taub-NUT geometry were analyzed using algebraic and analytical methods 
O HB1 E] . One of the purposes of this paper is to prove that the approach 
can be developed in order to include into consideration algebroid Taub-NUT 
backgrounds, in the context of the generalization of gauge-invariant theories 
fTE\ ITKj of the Dirac field. In a more explicit form, in the present work, we 
develop an algebroid 50(4, 1) gauge like theory of the Dirac fermions con- 
sidered for spacetimes with generic off-diagonal metrics, for instance, in the 
external field of the Kaluza-Klein monopole ^3] which is deformed to Lie al- 
gebroid confiugrations. Our aim is also to emphasize some new features of the 
Einstein theory in higher dimension spacetime when the locally anisotropic 
properties, induced by anholonomic constraints and extra dimension gravity, 
are characterized by Clifford algebroid symmetries. We construct new classes 
of exact solutions of the Einstein-Dirac equations defining 3D soliton-spinor 
configurations propagating self-consistently in a 5D Lie algebroid Taub NUT 
spacetime and analyze certain physical properties of such geometries. 

We note that in this paper the 5D spacetime is modelled as a direct time 
extension of a 4D Riemannian space provided with a corresponding spinor 
structure, i. e. our spinor constructions are not defined by some Clifford al- 
gebra associated to a 5D bilinear form but, for simplicity, they are considered 
to be extended from a spinor geometry defined for a 4D Riemannian space. 

We start in Section 2 with an introduction into the theory of (gravita- 
tional) nonholonomic Lie algebroids. In Section 3 we study 5D metrics char- 
acterized by Lie algebroid symmetries and nonholonomic constraints. Two 
examples of such exact solutions are analyzed. We also point the possibility 
to model corrections to the Newton low by extra-dimension polarizations 
and nonholnomic frames which is very similar to warped geometries. In Sec- 
tion 4, we consider the Dirac equations on gravitational algebroids. Then, 
we generate new solutions of the 5D Einstein - Dirac equations constructed 
as generalizations of algebroid Taub NUT vacuum metrics to configurations 
with Dirac spinor energy-momentum source, i. e. to Clifford algebroids. Fi- 
nally, in Section 5, we discuss and conclude the work. The Appendix contains 
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some ansatz formulas for Einstein equations with nonholonomic variables and 
their solutions. 

2 Spacetimes with Lie N— Algebroid Symme- 
try 

Let us consider a 5D spacetime V with a conventional splitting of dimension, 
dirriV = n + m = 5, where n > 2, which is a (pseudo) Riemannian mani- 
fold, or a more general one with nontrivial torsion (i.e. a Riemann-Cartan 
manifold) of necessary smooth class. The splitting can be defined in global, 
coordinate free form, as a locally non-integrable (nonholonomic) distribution 
globalized for every point u G V resulting in a Whitney type sum 

TV=hV®vV (1) 

stating a decomposition of the tangent bundle TV into certain conventional 
horizontal (h) and vertical (v) subspaces. We call a such manifold V to 
be N-anholonomic being provided with a nonlinear connection (in brief, N- 
connection) structure N defined by (JTJ), see details in Refs. [HI2|- 2 

We state the typical notations for abstract (coordinate) indices and geo- 
metrical objects defined with respect to an arbitrary local basis (system of 
reference, vielbein, or funfbein for 5D spacetimes) e a = (ei,v a ) on V. The 
small Greek indices a, j3, 7, ... run values 1, 2, . . . ,n + m and i,j, k, ... and 
a, b, c, ... respectively label the geometrical objects on the base and typical 
"fiber" and run, correspondingly, the values 1,2, ...,n and 1,2, ...,m. The 
dual base is denoted by e a = (e l , v a ). The local coordinates of a point u G V 
are written u =(x, u), or u a = (x l , u a ), where x l are local /i-coordinates, with 
respect to ei, and u a are local ^-coordinates with respect to the basis v a . 

A N-connection N is also given by its coefficients, 

N = N%u)dx i ®t^ = ® v b , 

where there are underlined the indices defined with respect to the local co- 
ordinate basis 

e« = da = d/du- = (ei_ = di = d/dx l , Vg, = da = d/du-) 

2 A manifold is called nonholonomic if it is provided with a nonintegrable distribution; 
in literature, one uses the equivalent term 'anholonomic'. In this papers, we consider 
a special case of nonholonomic manifolds when the anholonomy is defined only by the 
N-connection structure. 
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and its dual 

= du & = (e l = dx\ e a = du 3 ). 

The class of linear connections is parametrized by linear dependencies on u-, 
i. e. Nf (x, u) = r^(x)uK 

A N-connection structure induces a system of preferred vielbeins on V : 
Let us consider a 'vielbein' transform 

e a = e a ( u )e a and e a = e^(u)e^ (2) 

given respectively by a nondegenerated matrix e /3 ~(u) and its inverse e a /3 (u). 
Such matrices respectively parametrize maps from a local coordinate frame e« 
and co-frame e— to any general frame e a = (e,, v a ) and co-frame e a = (e\ v a ). 
If we consider a subclass of matrix transforms (0) linearly depending on 
N^x,^, with the coefficients 



e « (u 







(3) 



and 



e\(u) -iV fc & (u)e fc 4 (u) 



(4) 



we generate N-adapted frames 



(5) 



and their dual coframes 



e Q = ( e \ v fe = v b + N b i dx i ), (6) 

for any t> b = e b ~db satisfying the condition v c \v b = 5 b . In a particular case, we 
can take Vb = <9&. The operators © and (|BJ) are the so-called " N-elongated" 
partial derivatives and differentials which define a N-adapted differential 
calculus on N-anholonomic manifolds. 3 

The Lie algebroids structure on N-anholonomic manifolds (equivalently, 
Lie N-algebroids) is defined by the corresponding sets of functions fP a (x,u) 

3 we use 'boldface' symbols in order to emphasize that the geometrical objects are 
defined in N-adapted form, with invariant h- and v-components, on a manifold provided 
with N-connection structure 
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and C[ g (x,u), see details in Refs. [TJ |2j. For such Lie N-algebroids, the 
structure relations satisfy the conditions 

K v b) = Pb(x,u) e i; (7) 
[v d ,v b ] = C f db (x,u)v f (8) 

and the structure equations of the Lie N-algebroid are written 

te(/? 6 )-fe(/? a ) = f? e c e ab , (9) 
E (fe(c{j+c{ 9 cL-ci> a Q^;) = o, 

cyclic(a,b,e) 

for Q^blj = eb 'b ee ' e e f e j( e b~ e e ~ e ^ /) "with the values e 6 ' 6 and eyr defined by the 
N-connection. The anchor is defined as a map p : V — > /iV and the Lie 
bracket structure C^ 6 is considered on the spaces of sections S'ec(fV). For 
trivial N-connections, we can put iV" = and obtain the usual Lie algebroid 
constructions with eip\ — > dip\ when the structure functions p\{x) and C^ b (x) 
do not depend on v-variables u a . We can say that a Lie N-algebroid geometry 
is modelled on a spacetime V provided with nontrivial N-connection N and 
Lie N-algebroid structures p\ and C e ah subjected to the conditions (J7J) - Q. 
Such spacetimes are generalizations of the (pseudo) Riemannian manifolds 
with Lie group symmetry. 

The curvature of a N-connection Q = —N v is introduced as the Nijenhuis 
tensor 

N V (K, Y) = [ vX, vY] + vv [X, Y] - v[ vX, Y] - v[X, vY] 

for any vector fields X and Y on V associated to the vertical projection "v" 
defined by this N-connection, i. e. 

ft =ln b e l A e j <g> v b 
2 J 

with the coefficients 

= eyJV fl = - e,iV« + JV*^ (JV» ) - iV> 6 (N%) . 

The vielbeins (0) satisfy certain nonholonomy (equivalently, anholonomy) 
relations 

K,e p ] = W^e, (10) 
with nontrivial anholonomy coefficients 

W? k = Q$ k {x,u), Wl = v e N\(x,u) and W b ae = C h ae {x, u) 



7 



reflecting the fact that the Lie algebroid is N-anholonomic. 
A metric g on V can be written in N-adapted form, 

g = g Q/3 (u) e a <g> = 9ij (u) e l ® e> + h cb (u) v c ® v\ (11) 

where ^ = g(e i ,e i ) and h cb = g(v c ,v b ) and e u = (ei,v b ) and = (e\ v b ) 
are, respectively, just the vielbeins (jHJ) and (jHJ)- We can define the anchored 
map for the " contravariant" v-part of 

/i cfc (u) v c <g> v c — h cb (u) p* c p£ e, <g> Bj (12) 

modelling a h-metric ^/i 1 - 7 = /i cfe (u) pj, p^. There are certain anchors p| 
when ^/i 1 - 7 = (yf^. 

A distinguished connection (d-connection) D = {r§ } on V is a linear 
connection conserving under parallelism the Whitney sum (0). This mean 
that a d-connection D may be represented by h- and v-components in 
the form = (L) k , L bk , B l jc , B b ^j , stated with respect to N-elongated 

frames (jBJ) and (JHJ), defining a N-adapted splitting into h- and v-covariant 
derivatives, D = hD + vD, where hD = (L, L) and vD = (B, B). 

A distinguished tensor (in brief, d-tensor; for instance, a d-metric (jllj) ) 
formalism and d-covariant differential and integral calculus can be elabo- 
rated 13 IUE1 for spaces provided with general N-connection, d-connection 
and d-metric structure by using the mentioned type of N-elongated oper- 
ators. The simplest way to perform a d-tensor covariant calculus is to use 
N-adapted differential forms with the coefficients defined with respect to (ffij) 
and ©, for instance, T^j = Y^e 1 . 

The torsion T a and curvature lZ a * are defined by standard formulas but 
for N-adapted differential forms: we have respectively 

T° = De a = de a + T« A (13) 

and 

see Refs. [TJE] for explicit formulas T a = {T^} and TZ a p = {RV} for the 
coefficients computed with respect to N-adapted frames © and (jBj). 

The Ricci d-tensor R ^ can be computed by contracting the correspond- 
ing indices 

R f3-y =F 

and scalar curvature is _ 

R = g^R PT 
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The Einstein equations are written in the form 

G 7 = R p 7 — -g /3 7 R = T a/3 (15) 

for a general source, T a p, of matter fields and possible extra dimension cor- 
rections. 

A Riemann-Cartan algebroid (in brief, RC-algebroid) is a Lie algebroid 
A = (V, [•, •] , p) associated to a N-anholonomic spacetime V provided with 
a N-connection N, symmetric metric g(u) and linear connection T(u) struc- 
tures resulting in a metric compatible covariant derivative D, when Dg = 0, 
but, in general, with non-vanishing torsion. In this work, we shall inves- 
tigate some classes of metrics g(u) and linear connections T(u) modelling 
RC-algebroids as exact solutions of the 5D Einstein-Dirac equations. 

On RC-algebroids, the Levi-Civita linear connection V = { V T^ 7 }, by 
definition, satisfying the metricity and zero torsion conditions, is not adapted 
to the global splitting (|TJ) and can not applied for elaborating N-adapted and 
algebroid constructions, see details and discussion in Refs. [UE]- Neverthe- 
less, there is a preferred canonical d-connection structure V constructed only 
from the metric and N-connection coefficients [gij, h ab , iVf ] and satisfying the 
metricity conditions Dg = and T l - k = and T a bc = 0. In explicit form, the 
h-v-components of the canonical d-connection r 7 a/3 = (L* fc ,L£ fc , B % - Cl B bc ), 
are given by formulas 



p 


= \f '[e*(&>) + ej(9kr) - 




fa 
n bk 


= <N a k ) + hi ac [e k {h bc ) 


- h dc v b {N d ) - h db v c (N^)] , 


B) c 






B bc 


= hi ad [v c {h u ) + v b {h cd ) 


- v d (h bc )} . 



The formulas (fTUJl . (JT3j) and (JTljl are defined on the nonholonomic space- 
time V and contain the partial derivative operator v c = d/du c . We can 
emphasize the Lie N-algebroid structure by working with "boldface" oper- 
ators v c — > v c = p^x, u)Gi (see formulas (|ZJ) and ©)• A such "anchor- 
ing" of formulas defines canonical maps for d-metrics, anholonomic frames, 
d-connections and d-torsions from V to Sec(t>V). By anchoring the N- 
elongated differential operators, we can define and compute (substituting 
v c by p^ c <&i into (|T6|) ) the canonical d-connection p T 1 a/3 on Sec(vV) stating 
a canonical map T 7 Q/3 — * p ~T J al3 . 
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3 Algebroid Taub NUT spaces 



The standard Kaluza-Klein monopole was constructed by embedding the 
Taub-NUT gravitational instanton into 5D theory, adding the time coordi- 
nate in a trivial way ^T] . There were investigated locally anisotropic variants 
of such solutions when anisotropies are modelled by effective polarizations 
of the induced magnetic field. The aim of this Section is to analyze nonholo- 
nomic deformations of the Taub-NUT solutions when the metrics possess 
Lie N-algebroid symmetry. 

3.1 Background metrics and deformations to gravita- 
tional algebroids 

The Taub NUT solution of the 5D vacuum Einstein equations {R a p = 0, for 
the Levi-Civita connection) is expressed by the line element 

,2 



ds {5D) = dt +ds {4D) , (17) 

\AD) 

where 



dsf m = -Q-\dr 2 + r 2 d9 2 + sin 2 6d<p 2 ) - Q{dx 4 + Aidx 1 ) 2 



Q 1 = 1 H -,Tn = const. (18) 

r 

The functions are static ones associated to the electromagnetic potential 

A r = 0, A g = 0, A v = 4m (1 - cos 9) 
resulting into "pure" magnetic field 

B = iotA = m ^j (19) 

of a Euclidean instanton; is the spherical coordinate's unity vector. 

The metric ()17j) defines a spacetime with global symmetry of the group 
G s = SO(3) <8> £^4(1) ® T t (l) since the line element is invariant under the 
global rotations of the Cartesian space coordinates and x 4 and t transla- 
tions of the Abelian groups ^4(1) and T t (l) respectively. We note that the 
£7,4(1) symmetry eliminates the so called NUT singularity if x 4 has the period 
Anm . The mentioned group symmetries can be deformed in Lie N-algebroid 
one by corresponding nonholonomic frame transforms. Let us consider this 
procedure in details: 

We introduce a new 5th coordinate, 

x ^ q = x ^ f ^(9,^(9^), (20) 
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when 

dq + 4m (l — cos 6)d6 = dx 4 + 4m (l — cos 6) dip. 
This holds, for instance, for 

di = n(0, ip)d(, - y 4 ) = ^d6 + ^dip, 

when 

— = 4m (l — cos 6>)/i and — — = — 4m (l — cos#)/x 
ad dip 



with 



/x = (1 — cosO) 2 exp [9 — ip] . 



The changing of coordinate ()20|) describes a reorientation of the 5th coordi- 
nate in a such way as we could have only one nonvanishing component of the 
electromagnetic potential 

A e = 4m (1 - cos6») . 
The next step, we consider an auxiliary 5D metric 4 



ds%D) = dt 2 + ds 2 (m , (21) 
dsf 4D) = -dr 2 - r 2 d9 2 - r 2 sin 2 6dip 2 - Q 2 (d( + A e d6) 2 , 

where 

This metric will be transform into some exact solutions after corresponding 
N-anholonomic transforms. 

Let us consider a 5D ansatz of type (jllj) 

Ss 2 = (dx 1 ) 2 + g 2 (x k ) (dx 2 ) 2 + g 3 (x k ) (dx 3 ) 2 
+h 4 (x k ,v)(e i ) 2 + h 5 (x k ,v)(e 5 ) 2 , 
e 4 = dy 4 + Wi(x h , v)dx i and e 5 = dy 5 + rii(x k , v)dx i (22) 

with the time like coordinate x 1 = t and "anisotropic" dependence on coor- 
dinate y 4 = v and running of indices like ... = 1, 2, 3 and a, b, ... = 4, 5. 
The set of coordinates x k = (x 2 , x 3 ) and y 5 can be any parametrization of 
the space line and the 5th dimension coordinates. This way, the complete 

4 it is not a solution of the Einstein equations, but its corresponding deformations will 
generate a class of exact solutions 
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set of local coordinates is stated in the form u = {u a = (t,x 2 ,x 3 ,y i ,y 5 )}. 
The N-connection coefficients are parametrized in the form Nf = Wi(x k ,v) 
and Nf = rii(x k ,v). We also write 

g t = q i (x k )7 li (x k ) and h a = q a {x k )r] a {x k , v) (23) 

where r\ a = (1,7ft, r] a ) are called 'polarization' functions. The values q a = 
(1, (ft, q a ) are just the coefficients of the metric (J2~T|) if 7] a — > 1 and u>i, — ► 0. 

Our aim is to find certain nontrivial values r] a and Wi, rti when the metric 
(|22jl defines a solution with Lie N-algebroid symmetry, i. e. there are satisfied 
some conditions of type (fl2~|). 

0«(u) = /i«*(u) £ (u) ft(u). 

Such anchor conditions, for effectively diagonalized metrics (with respect to 
N-adapted frames), must be satisfied both for rj a — 1 and nontrivial values 
of T] a , i. .e. 

<? = {f? 4 ) 2 + h 5 (/? 5 ) 2 , forr/^1; (24) 
q* = g 4 (tiY + q 5 {f?,)\ for ry a = 1, 

where = l/gi,h a = l/h a ,r] a = l/r) a and g° = 

By straightforward computations (see explicit formulas and details in 
Refs. P3 12 IU E]) we can check the nontrivial components of the Einstein 
tensor G°jg for the canonical d-connection r 7 a/3 (fTBj) satisfy the conditions 

G{ = -(Rj + Ri), Gl = G\ = -R\{x\ x\ v),G 4 4 = Gl = -R 2 2 (x\ x 3 ). 

This mean that the Einstein equations (JT5|l for the ansatz f)22j) are compatible 
for nonvanishing sources and if and only if the nontrivial components of the 
source, with respect to the frames © and (jBj) . are any functions of type 

fj = f 3 = T 2 (x 2 , x 3 , v), T4 = Tg = T 4 (x 2 , x 3 ) and f\ = T 2 + T 4 . (25) 

Parametrizations of sources in the form (|23|) can be satisfied for quite general 
distributions of matter, torsion and dilatonic fields in string gravity or other 
gravity models. In this paper, we shall consider that there are given certain 
values T 2 and T 4 which vanish in the vacuum cases or can be induced by 
certain packages of spinor waves. 

A very general class of exact solutions of the Einstein equations (JXEJ), 
with nontrivial sources of type fl25(l parametrized by the metric ansatz f!22D . 
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see Appendix (for simplicity, those formulas derived for the conditions (|59|T 
is described by off-diagonal metrics of type 

Ss 2 = (dx 1 ) 2 - g~ k (/) (dx k f - 

h 2 0] (*<) [/' (x\v)] 2 |, x (*>) I (e 4 ) 2 - [/ (x\v) - / (xO] 2 (e 5 ) 



2 



e 



! dv + w k (x\ v) dx k , e 5 = dv 5 + n k {x\ v) dx k , (26) 



where the coefficients yx l j are constrained to be a solution of the 2D 
equation 

.. _ gjgi _ igr + <> _ ?A _ (M_ = 2gmTi{x ^ x3) (27) 

2^2 2^3 2# 3 2^2 

for a given source T 4 (^x l ^j where, for instance, g° = dg-^/dx 2 and g' 2 = 

dg 2 /dx 3 and /* = df /dv. It is always possible to find solutions of this equa- 
tion, defining 2D Riemannian metrics, which are conformally flat, at least 
in non-explicit form. Hereafter we shall consider that g k (^x l ^j are certain 

defined functions. The rest of functions from ()26j) can be computed in the 
form: 

, r (x\v) = 1 - 1 J T 2 (x % ,v)-^[f (x\v) - f {x l )] 3 dv; (28) 
the N-connection coefficients Nf = Wi(x k ,v) and Nf = rii(x k ,v) are 

C T {x\v) [ZJ) 

and 

r d 

n k = n k[1] (x l ) + n m (x l ) / ^ T (x\ v) — [f (V , v) - f (x*)} 3 dv. (30) 

The set of functions (j2*Tj) - (j3TH) defines a class of exact solution of the 
5D Einstein equations depending on arbitrary nontrivial functions f(x\v) 
(with /* 7^ 0), hl(x l ), 54[ ] (x % ) , n k [i] (x*) and n k [ 2 ] (x l ) , and sources T 2 (x k ,v) 
and T 4 (^x l ^j which have to be defined from certain boundary conditions and 
physical considerations. 5 It is not difficult to see that the metric ()26|) defines 



5 Any metric l)26|l with h\ ^ and /15 =^ has the property to be generated by a func- 
tion of four variables / (x l ,v) with emphasized dependence on the anisotropic coordinate 

v, because /* = d v f ^ and by arbitrary sources T2(x^,v), T4 ■ The rest of arbi- 
trary functions not depending on v have been obtained in result of integration of partial 
differential equations. 
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a Lie N-anholonomc algebroid. This follows from a parametrization of type 

9t = QiVh -h[ 0] {x l ) [f*] 2 \<;r\ = qtfk, - [f - /of = ?5?75, (31) 

which allows to define the polarization functions r] a and compute the non- 
trivial anchor coefficients (|24j) . 

The Lie N-algebroid structure is finally defined for such classes of metrics 
if the structure functions C^ b (x, u) are chosen to satisfy the algebraic relations 
with given values for (|2*3j) an d defined N-elongated operators e^. In 
result, the second equation in © will be satisfied as a consequence of the 
first one. This restricts the classes of possible v-frames, v\, = e b ~(x,u)d/du-, 
where e b ~(x,u) must solve the algebraic equations (jHJ). We conclude, that 
the Lie N-algebroid symmetry imposes certain algebraic constraints on the 
coefficients of vielbein deformations generating the off-diagonal solutions. 

The sourceless case with vanishing T 2 and T 4 can be distinguished in 
the form: Any off-diagonal metric ([26)) with h^x 1 ) = = const, Wi = 
and rik computed as in (}30j) but for $y — 1 defines a vacuum solution of 
5D Einstein equations for the canonical d-connection pfi)l . By imposing 
additional constraints on arbitrary functions from Nf = rii and Nf — Wi, we 
can select just those off-diagonal gravitational configurations when the Levi- 
Civita connection and the canonical d-connections are related to the same 
class of solutions of the vacuum Einstein equations, see details in Ref. [H El 
ElEj. With the aim to analyze the gravitational algebroids in general form, 
in this paper, we shall consider nontrivial torsion configurations induced by 
nonholonomic frames. 

3.2 Two examples of Taub NUT algebroids 

We outline two classes of exact solutions of 5D vacuum Einstein equations 
on Lie N-algebroids which, in Section |4~21 will be extended to configurations 
with spinor matter field source. 

3.2.1 Static gravitational algebroids with angular polarization 

A stationary ansatz for of type (f2*Uj) with explicit dependence on the "aniso- 
tropic" angular coordinate v = if is taken in the form 

dt 2 — 5s 2 AD ^ , 

-dr 2 - r 2 d9 2 - r 2 sin 2 9 r] 4 (r, 6, <p)5<p 2 - Q 2 (r) rj 5 (r, 6, <^)ck 2 , 
dip + u>i(r, 9, ip)dt + W2(r, 9, ip)dr + ws(r, 9, ip)d9, 
dq + ni(r, 9, <p)dt + n 2 {r 1 9, ip)dr + n 3 (r, 9, <p)d6, 



5s 2 = 

^s 2 

0S {iD) — 

5ip = 

<k = 
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where g 2 = — 1 and q 3 = — r 2 solves the equation (J27j) for T 4 = and 
7?4,5 ( r ; 0> V 9 ) define classes of exact vacuum solutions with T 2 = 0. For sta- 
tionary metrics, one does not consider dependencies of the metric coefficients 
on the time coordinate. 

The formulas PKJ) and (J2HI), for Y 2 = 0, allow respectively to take $y — 1 
and W{ = 0. We also consider ni = and n 3 = which can be obtained from 
(JHUj) by stating zero the corresponding integration functions, i.e. niji^j = 
and ri3[i,2] = 0. The formulas (|3*T|) give 

g 2 = -1, r? 2 = 1; <?3 = -^ 2 , ??3 = 1; ?4 = sin 2 0, g 5 = -Q 2 (r) (32) 

where 774 (r, 6 1 , <p) and 775 (r, 97) satisfy the conditions 

frfo] ( r ) LTfr = r2 sin2 # ^(r, 9, tp), [f(r, tp) - f (r)} 2 = Q 2 (r)r] 5 (r, tp). 

(33) 

Such conditions are satisfied for any set of functions f(r, <p), fo(r), fy ] (r, 9) , 
774(7", 9, ip) and 775(7", tp) having limits 77 4)5 — ► 1 for <p — > 0. 

In the locally isotropic limit of the solution for n 2 (r, 0, <p) , when <£> — > 0, 
we obtain the particular magnetic configuration contained in the metric (|2*T|) 
if we impose the condition that 

f d 

n 2 [o] + n 2 [ii hm / — [/ (r, 0, <p) - / (r, 0)]~ 3 dv9 = A 9 = 4m (1 - cos0) , 
tp^O j otp 

which defines only one function from two unknown values n 2 [o] (r, 9) and 
n 2[i] { r i 9) ■ This has a corresponding physical explanation. From the usual 
Kaluza-Klein procedure, we induce the 4D gravitational field (metric) and 
4D electromagnetic field (potentials Ai), which satisfy the Maxwell equa- 
tions in 4D pseudo-Riemannian spacetime. For the case of spherical, locally 
isotropic, symmetries the Maxwell equations can be written for vacuum mag- 
netic fields without any polarizations. When we introduce into consideration 
anholonomic constraints and non-spherical symmetries, the effective mag- 
netic field could be effectively polarized by higher dimension gravitational 
field or vacuum nonlinear gravitational interactions. In the simplest case, we 
can put n 2 [oj (r, 9) = and n 2 [i] (r, 9) = 4mo (1 — cos0) and write 

n 2 (r, 9, tp) = Am (r, 9, tp) (1 — cos 9) 

where the gravitationally anisotropically polarized mass is defined 

d 



m(r,9,<p)= I ^[/(r,^)-/oM)]- 3 ^, 
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with 

lim m (r, 9, (p) = m . 

95 — >o 

This class of 5D vacuum gravitational static algebroid solutions can be 
represents in the form 

Ss 2 = dt 2 -5s 2 (4D) , (34) 
Ss 2 (4D) = ~dr 2 - r 2 d6 2 - r 2 sin 2 9 77 4 (r, 6, <p)d<p 2 - Q 2 {r) T] 5 (r, 6, y?)ck 2 , 
5q = dq + 4m (r, 9, tp) (1 — cos 9) dr. 

Considering r\\ = 1 and q\ = 1 and introducing the gravitational polariza- 
tions r] a for the data and (jSSJ), we can compute the anchor coefficients by 
solving the algebraic equations (|2*4^1 and define the Lie N-algebroid structure 
with any C£ c satisfying the conditions Q. Such d-metrics are similar to the 
Taub NUT vacuum metric (fTTj) and its 4D conformally transformed partner 
(j21J) with that difference that the coefficients are polarized by nonholonomic 
constraints. Their Lie algebroid symmetry is a nonholonomic defibrinations 
of the global symmetry defined by the group G s = SO(3) <8> U^il) <£> T t (l). 
We can treat (|H4^1 as a static algebroid Taub NUT solution. 



3.2.2 Solutions with extra— dimension induced polarization 

Another class of solutions is constructed if we consider a d-metric of the type 
()26|) with explicit dependence on the "anisotropic" angular coordinate v = q 
is taken in the form 

dt 2 — 5S(4D), 

-dr 2 - r 2 d9 2 - Q 2 {r) T] 4 (t, r, 9, ^)<k 2 - r 2 sin 2 9 T] 5 (t, r, 9, <;)5<p 2 , 
dq + Wi(t, r, 9, q)dt + w 2 (t, r, 9, q)dr + w 3 (t, r, 9, q)d9, 
dip + m(t, r, 9, q)dt + n 2 (t, r, 9, q)dr + n 3 (t, r, 9, q)d9. 

The equation (|2Tj) is satisfied for T 4 = 0, but 774,5 (t, r, 9, are defined for a 
class of exact solutions with nonzero source, T 2 (r, 9) 7^ 0, in equation (jo^j) 
from Appendix. 

For simplicity, we can consider stationary solutions when the functions 
7/4,5 and Wi,rii do not depend on time variable t and consider functions of 
type f(r,9,q),f (r,9) and any integrable T 2 (r, 9). The solutions (|2~%j). (}2T?j) 
and ()30|) are respectively written in the form 

?r (r,9^) = l-^J T 2 (r, 0, ? )^[/ (r, 0, - /„ (r, fl)] 3 dq 



5s 2 = 

<5 S (4L>) — 

<k = 

5p = 
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the N-connection coefficients Nf = Wi(x k ,v) and Nf = rii(x k ,v) are 

<9<; T (r, 6, ?) / (0, ^) 
wi = 0, w 2 = — 



w 3 



dr \ 
d^r (r, 0, / <9^ T (0, \ _i 



30 V 
and 

^ T (r, 0, f ) — [/ (r, 0, f ) - / (r, 0)] _3 tfc. 

We have to impose the limits 

w 2 — > and u> 3 — >■ A e = 4m (l — cos0) for c — >■ 0, 

this can be obtained by a corresponding class of functions T 2 (0, and 
/ (0, <;) , f (0) such that 

lim[/ (r, 0, s ) - f Q (r, 0)] 2 = r 2 sin 2 for lim[/* (r, 0, q)} 2 = Q 2 {r) 

and nfc — > for <j — > in order to get (|7T|) . 
Expressing 

l^| = ^(r,0)[/>,0,<)] 2 |<; T (r, and |/i 5 | = [/(r, 0,?) - f (r, 0)] 2 , 
^4 = [r(r,0^)] 2 andr / 5 = [l-/(r,0,O//o(r,0)] 2 (35) 

and parametrizing the integration functions to have 

= 1 - y 2 J T 2 (r, 0, (r, 0, ,) - / (r, 0)] 3 , (36) 

«;! = o,T 02 = i/*|- 1 |:i/-/oi,ws = irr l ^i/-/oi, 

n k = n k [x] (r,0) + n k[2] (r, 0) (/ - / )~ 2 , 
with the limits / — > 0, /* — >■ 1 and 

irr 1 ^l/-/o|-4m o (l-cos0) 

for <j — > and hl(r,8) = Q 2 (r) = — g 4 and (/ ) 2 = r 2 sin 2 = — g 5 . We 
have nfe — > with / — > for ^ — > if n^fi] (r, 0) = —n k \^ (r, 0) . We put 
ni[i] (r, 0) , nip] (r, 0) = in order to get static solutions. 
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The data for a such solution are concluded for the metric 



Ss 2 = dt 2 -5s 2 m , (37) 



Ss 2 4D) = -dr 2 -r 2 d9 2 



-Q 2 (r) [/*(r,0,O]VM,Olfc a - [rfrin0-/(r,0 J ?)] S V, 
ft = rf^+l/T 1 ^\f-Mdr+^\f-f \d9 , 
5ip = dip + [1 - (/ - / r 2 ] (r, 0)dr + n 2[2] (r, 9) 



generated by functions f(r,9,q),fo(r,9) and n 2 [i](r, 9), n 2 [ 2 ](r, 9) satisfying 
the above stated conditions. In this case we have r\\ — 1 and qi — 1 : intro- 
ducing the gravitational polarizations r/a for the data (|3~T|) and (pT7jl . we can 
compute the anchor coefficients by solving the algebraic equations ({231 and 
define the Lie N-algebroid structure with any C£ c satisfying the conditions 
• We note that this metric can not be transformed into a vacuum because 
the (jHHJ) with W3 7^ are possible for 7^ 1 and T 2 7^ 0. Such type of 
gravitational algebroid configurations can be generated from the Taub NUT 
solution by nonholonomic transforms only by nontrivial matter sources or 
by any matter like corrections from extra dimension, for instance, in string 
theory. The polarization functions depend explicitly on extra coordinate. 

Finally, we analyze the possibility to generate warped (on the extra 
dimension) gravitational algebroid configurations. In the original Randal- 
Sundrum scenaria the Newtonian gravitational potential takes the form 

V[t) = G N 



r 2 k 2 

for two interacting point masses mi and m 2 , where Gn is the 4D gravitational 
constant and k is the "warping" factor. In the papers jUl Ej we concluded 
that anholonomic coordinates can give rise to " anisotropic" deviations of the 
Newton potential, like 



V{r) = G 



171x1712 



N- 



,-2k y \y\ 

r 2 k 2 



where y can be a space like or extra dimension coordinate. We can consider 
such deviations in the Taub NUT algebroids. For instance, we consider 
instead of Q _1 ((IHI) a function 

Q- 1 = 1 + (l + ^J ,m ,k = const. (38) 
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It is possible to substitute Q — > Q and/or to include a factor of type (1+ 
e~ 2k v\y\/r 2 k 2 ) into 775 for the class of solutions (|34p. In result, one gener- 
ates warped static gravitational algebroids from a vacuum solution but with 
nonholonomic (off-diagonal) metric terms which play the role of source from 
the "isotropic" brane constructions. More similarities with the former brane 
constructions may be considered for the class of metrics (|37|) where T 2 can 
be approximated as a constant tension (source of anisotropy), in general de- 
pending on extra dimension, for the second class of graviational algebroids. 



4 Einstein— Dirac Algebroids 

The geometry of spinors on Clifford algebroids is elaborated in Ref . |T] . Such 
spinors can be included in the left or right minimal ideals of Clifford alge- 
bras generalized to Clifford spaces (C-spaces, elaborated in details in Refs. 
[10j), in our case provided with additonal N-connections and/or algebroid 
structure. In this Section, we consider Dirac spinors defined with respect to 
N-anholonomic frames on 5D gravitational algebroid spaces. 

4.1 Dirac equations on gravitational algebroids 

For a d-metric ()22|) with coefficients 

9 a p(u) = (9ij(u),h ab (u)) = (l,^(w),/i a (w)), 

where i = 1,2; i = 0,1, 2; a = 3,4, defined with respect to an N-adapted 
basis we can easily define the funfbein (pentad) fields 

e E = = {e tL = e%, e« = e^dj, (39) 

ei± = e ^ = {e i = e\d\ = e%5 a } 

satisfying the conditions 

gij = diag[l, —1 — 1] and hab = diag[—l, —1]. 
The d-metric (|2*2*|) is effectively diagonal. This allows to write 

4 = V\9T\SI and ej = 
where 5f and 8% are Kronecker's symbols. 
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The Dirac spinor fields on nonholonomically deformed Taub NUT spaces 

y(u) = iV a (u)] = [^ (u),xt(u)}, 

where / = 0, 1, are defined with respect to the 4D Euclidean tangent subspace 
belonging the tangent space to the 5D N-anholonomic manifold V. The 4x4 
dimensional gamma matrices 7- = [7- , 7- , 7- , 7- ] are defined in the usual 
way, in order to satisfy the relation 

{^',^}=2g^, (40) 

where ^7- 7- | is a symmetric commutator, g--^ = (—1,-1,— 1,-1), which 

generates a Clifford algebra distinguished on two holonomic and two anholo- 
nomic directions. In order to extend the ()40|) relations for unprimed indices 
a, (3... we conventionally complete the set of primed gamma matrices with a 
matrix 7-, i. .e. write 7- = [7-, 7-, 7-, 7-, 7-] when 

{7^, 7^} = 

The coefficients of N-anholonomic gamma matrices can be computed with 
respect to anholonomic bases (JSJ) by using respectively the funfbein coeffi- 
cients 

We can also define an equivalent covariant derivation of the Dirac spinor 
field, V a $, by using pentad decompositions of the d-metric (j22J), 

(41) 

where there are introduced N-elongated partial derivatives and the coeffi- 
cients 

(Sa^y (u) = (D 7 e^) e^e^ 

are transformed into rotation Ricci d-coeflicients Sap-y which together with 
the d-covariant derivative D 7 are defined by anholonomic pentads and an- 
holonomic transforms of the Christoffel symbols. In the canonical case, we 
should take the operator of canonical d-connection D 7 with coefficients (fTo^) . 

For diagonal d-metrics, the funfbein coefficients can be taken in their turn 
in diagonal form and the corresponding gamma matrix 7" (u) for anisotropic 
curved spaces are proportional to the usual gamma matrix in flat spaces 7X 
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The Dirac equations on Clifford algebroids [B] are written in the simplest 
form with respect to anholonomic frames, 

(ry a (u)vl-^ = 0, (42) 

where \x is the mass constant of the Dirac field. The Dirac equations are the 
Euler-Lagrange equations for the Lagrangian 

C (1 ^ (u) = («) 7° («) V> (u) (43) 

-(V> + (u))T (u) (u)] - fiV + (u) tf («)}, 

where by (w) we denote the complex conjugation and transposition of 
the column \l/ (u) . Varying Z^ 1 / 2 ) on d-metric (|4Hjl we obtain the symmetric 
energy-momentum d-tensor 

?«p(u) = ^[vI/ + (n)7 Q (n)V^M/(n) + M/ + ( M )^(n)^H 

-(V>+ (u))% (u) * (u) - (V>+ (u))7„ («) * («)]. (44) 

By straightforward calculations we can verify that because the conditions 
D 7 e° = are satisfied the Ricci rotation coefficients vanishes, 

S«/3 7 (u) = and V Q # = 

and the N-anholonomic Dirac equations (}42|) transform into 

(«7 a («)e Q -^ = 0. (45) 

Further simplifications are possible for Dirac fields depending only on 
coordinates (t,x 2 = r,x 3 = 8), i. e. \& = \I/(x fe ) when the equation (|43|) 
transforms into 

{i^d t + %^^=d 2 + i^— l =dz - n)V = 0. 

The equation (j43j) simplifies substantially in ^-coordinates 

(t,C 2 = C 2 (r,^),C 3 = C 3 (r^)), 
defined as to be satisfied the conditions 

<9 2 and = (9 3 (46) 



1^2 1 a S V 03 
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We can consider a more simple equation 







d 



_d_ 



/i)*(t,C 2 ,C 3 )=o. 



(47) 



The equation (jlTf) describes the wave function of a Dirac particle of mass \i 
propagating in a three dimensional Minkowski flat plane which is imbedded 
as an N-adapted distribution into a 5D Lie N-algebroid. 

The solution ^ = ty(t, ( 2 , £ 3 ) of ()47|) is searched in the form 



for positive energy; 

v&(-)(C) = exp [zfaf + *2C 2 + fc3C 3 )k 1 (*0 
for negative energy, 



with the condition that k\ is identified with the positive energy and ^(k) 
and x 1 ^) are constant bispinors. To satisfy the Klein-Gordon equation we 
must have 

k 2 = (h) 2 - (k 2 ) 2 - %f 
The Dirac equations implies 



(a l ki - [i)(p l (k] and (a l ki + A*)x 1 (fc), 

where a % [i = 1,2,3) are Pauli matrices corresponding to a realization of 
gamma matrices as to a form of splitting to usual Pauli equations for the 
bispinors ^(k) and x l {k>)- 

In the rest frame for the horizontal plane parametrized by coordinates 
( = {t, C 2 , C 3 } there are four independent solutions of the Dirac equations, 



¥>(i)0,o) 



X(i)0,o) 



/ i \ 




V o y 



1 

V o y 



¥>(2)(M) 



X(2)(^0) 



1 



v o y 




V 1 / 



We consider wave packets of type (for simplicity, we can use only super- 
positions of positive energy solutions) 



d 3 p 



27T 3 ^2 + ( fc 2)2 



^ 6(p,N)^ N WexpHA;,Ci (48) 



>]=1,2,3 
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when the coefficients b(p, [a]) define a current (the group velocity) 




27r3 vV 2 + (k 2 ) 



d 3 p fj, 



f |6(p, N)l 2 



with \p 2 \ ~ and the energy-momentum d-tensor (|44|) has nontrivial coef- 
ficients 



T\ = 2T(C 2 , C 3 ) = fc^+Ta*, T 2 = T 3 = (49) 



where the holonomic coordinates can be reexpressed C — Ci xt )- We must 
take two or more waves in the packet and choose such coefficients b(p, [a]), 
satisfying corresponding algebraic equations, in order to get from (}4*9*j) the 



required by the conditions ([44)1 . 

Finally, in this Section, we note that the ansatz for the 5D metric (}2*2*|) and 
4D spinor fields depending on 3D h-coordinates £ = {t, ( 2 , ( 3 } reduce the 
Dirac equations to the usual ones projected on a flat 3D spacetime. This con- 
figuration is N-adapted, because all coefficients are computed with respect 
to N-adapted frames. The spinor sourse T(x 2 ,x 3 ) induces a corresponding 
Clifford algebroid configuration. 



In this subsection, we construct two new classes of solutions of the Einstein- 
Dirac fields generalizing the Taub NUT metrics defined by data ()34)1 and 
(J3*7j) to be respective solutions of the Einstein equations, see in the 

Appendix, with a nonvanishing diagonal energy momentum d-tensor 



for a Dirac wave packet satisfying the conditions ([49)1 and (|50|). 

4.2.1 Clifford algebroids with angular polarizations 

The vacuum d-metric (|34|) was constructed by taking qr = 1 for T2 = 
0(J2EJ) PSj) . For a nontrivial Dirac spinor source T 2 = T 2 = T(r,9) (f5U|). 
we compute a nontrivial matter polarization <^x defined in general form by 
introducing this source in formula (|28|). This results in nonzero values of Wj, 
defined by c^r in formulas (|2l?jl and modified rii because is also present 
in the formulas (|3())1. The integration functions in Wj and rij can be any 



equalities 



T 2 = T 3 = T(C 2 ,C 3 ) = T(x 2 ,x 3 ), 



(50) 



4.2 Algebroid Taub NUT 



Dirac Fields 



T2 = {2T(r,0),T(r,0),T(r,0),O,O} 
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ones subjected to the conditions that W\ — 0,n\ = 0, W7 2>3 — ► 0,n 3 — > but 
n 2 — > 4m (r, 9, ip) (1 — cos 9) for T — * 0. For a such source, the d-metric (J54*)l 
transforms into the form 

Ss 2 = dt 2 -5s 2 4D) , (51) 
5s 2 m = -dr 2 -r 2 d6 2 

-r 2 sin 2 9 r] 4 (r,9,p)\q T (r,9,(p)\Sp 2 - Q 2 {r) r]r (r,9,p)Sq 2 , 
5ip = dip + w 2 (r,9,{p)dr + w 3 (r,9,(p)d9, 
5q = dq + n 2 (r,9,(f)dr + n 3 (r,9,p)d9, 

where 

^ T (r,9,<p) = 1 - ^ J T(r, 9)^-[f (r, 9, <p) - f (r, 9)fdp. 

The nontrivial N-connection coefficients are computed 

d r sr (r, 9, ip) d e qr (r, g, y?) 

J <9^ T (r, 0, v?) ' 3 <%,9r (r, 9, tp) ' 

for d r = d/dr, d e = d/dO, d 9 = d/d<p, and 

n 2 ,3 (r, 9, ip) = n 2 , 3 [i] (r, 9) + n 2)3 [2] (r, 0) x 

y ^ (r, 0, |- [/ (r, 0, - / (r, 0)]" 3 ^. 

The data ([53]) are also modified because of q^ 

h [o] ( r ) [/*(**, ¥>)] 2 kx(r, 0,<p)\ = r 2 sm 2 9 rn(r,9,(p). 

This results in a modifications of the Lie algebroid anchor structure functions 
because (^—modified 774 and /14 (see formula (J2BJ)) change the solution of the 
algebric equations for ()24j1 . So, the anchor structure is also deformed if a 
vacuum gravitational algebroid is deformed to a such type of Einstein-Dirac 
algebroid. 

4.2.2 Clifford algebroids with extra dimension polarizations 

The solution (|37|) was constructed for a general matter source T 2 (r, 9,q). 
The spinor source (|5U|) can be considered a particular case when the matter 
energy-momentum tensor does not depend on extra dimension coordinate q, 
i.e. T 2 = T(r, 9), which results in a particular type of polarization qy(r, 9, q) 
computed just for a such T 2 (r, 9), see data (jBEJ). So, the d-metric (fHTj) 
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describes the solutions of the Einstein-Dirac equations as particular cases 
characterized by a proper anchor configuration because this type of Sr( r , Q, 
is related to (J35|) and (|2*1|) defining the algebroid configuration. This type of 
Clifford algebroids are also static but with coefficients depending explicitly 
on extra dimension coordinate <j. Such d-metrics are determined by the 
Dirac spinor field and do not have limits to vacuum configurations. 

Finally, we emphasize that all types of solutions considered in this work 
can be generalized to stationary configurations by introducing certain coef- 
ficients w\dt and /or n\dt in the off-diagonal part. Such terms have to be 
computed by corresponding formulas (|fi2*|) and with nontrivial integra- 
tion functions depending on h-coordinates. 

5 Discussion and Conclusions 

We have constructed a new class of nonholonmically deformed Taub NUT 
spacetimes possessing Lie algebroid symmetry. Such static gravitational al- 
gebroid configurations were generalized for nontrivial sources of Dirac spinor 
fields, i. e. the metrics were extended to define exact solutions of the 
Einstein-Dirac equations. They consist explicit examples of nonholonomic 
manifolds provided with distributions defining nonlinear connections when 
the Lie algebra symmetry was deformed to a Lie algebroid or Clifford alge- 
broid symmetry. 

There were distinguished two classes of 5D algebroid spacetimes: The 
first one is stated by solutions with angular local anisotropy and can ex- 
tended from vacuum configurations to nonlinear gravitational spinor interac- 
tions. The second one describes nonholonomic gravitational configurations 
induced by spinor sources and depends explicitly on extra dimension coor- 
dinate. Such metrics do not have vacuum limits if the nonlinear connection 
and Lie algebroid structure functions are not trivial. All constructed classes 
of solutions have smooth limits to the usual Taub NUT metric, can be ex- 
tended to stationary configurations and further deformed to other algebroid 
or non-algebroid symmetries, for instance, to ellipsoid or toroidal config- 
urations, deformed to wormholes, with, or not, warped factors and/or 3D 
solitonic gravitational waves like was proven in Refs. jE] . 

Let us now conclude the properties of such Einstein-Dirac algebroids. 

1. They are described by generic off-diagonal metrics and related nonholo- 
nomic vielbeins with associated nonlinear connection structure. Such 
Einstein, Einstein-Cartan and Einstein-Dirac spacetimes are with po- 
larization of constants and metric and connection coefficients. In gen- 
eral, there are nontrivial torsion coefficients, induced by frame non- 
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holonomy. By imposing certain constrains the constructions can be 
transformed into 4D (pseudo) Riemannian configurations. 

2. The metrics and connections for the solutions possessing algebroid sym- 
metries depends on certain classes of integration functions. This is a 
typical property for the algebroid approaches to strings in gravity |H|: 
It follows from general properties of the systems of partial nonlinear 
equations to which the Einstein equations are reduced for generic off- 
diagonal ansatz. In the 'simplest' case of diagonal ansatz depending on 
one variable (like the Schwarzshild metric) the gravitational field equa- 
tions transform into a nonlinear second order differential equation. Its 
general solution depends on integration constants which are physically 
defined from certain symmetry and boundary conditions in order to 
get in the limit just the Newton potential for a component of metric. 
By applying the anholonmic frame method we construct more general 
classes of solutions derived for partial differential equations. Such so- 
lutions contain integration functions depending on two, three or four 
variables and the corresponding spacetimes possess symmetries charac- 
terized not only by structure constants but also by structure functions 
and nonholonomic distributions. 

3. We can restrict a set of vacuum or nonvacuum off-diagonal metrics by 
imposing additional physical conditions like the Lie algebroid symme- 
try and fixing certain limits to well known solutions (for instance, to the 
Taub NUT metric). Nevertheless, even in such certain depen- 
dence on some integrability functions is preserved. It can be eliminated 
only by fixing an exact system of reference with a prescribed type of 
nonlinear connection and structure anchor and Lie algebra type struc- 
ture functions. This fixes the type of gravitational polarizations in the 
vicinity of some points on a nonholnomic manifold. But the general 
conclusion is that: when we deal with gravitational algebroids, we op- 
erate with classes of metrics and connections and corresponding classes 
of symmetries, i. e. with sets of spacetimes. 

4. As explicit examples, we defined and analyzed some Clifford algebroid 
configurations describing packages of 3D Dirac waves self-consistently 
propagating in 5D nonholonomic spacetimes. They are distinguished 
with respect to the so-called N-adapted frames where the symmetry 
properties can be defined in the simplest way. With respect to local 
coordinate frames, the nonlinear gravitational-spinor interactions mix 
into generic off-diagonal metric and vielbein configurations depending 
on four coordinates. 
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5. In a more general context, the method can be applied for definition 
of algebroid solutions with additional noncommutative symmetries, 
supersymmetric and/or complex variables, quantum deformations in 
string gravity and for generalized Finsler spaces like we elaborated for 
'non-algebroid' but also nonholonomic spaces |U El El HE] • 
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A Einstein Equations and Anholonomic Va- 
riables 

We outline some necessary formulas for the nontrivial components of the 
Einstein equations (fT3|) for the d-metric (}2*2*|) and canonical d-connection 
(fTBj). see references [U 12 IH E] for details on proofs and related references. 
One obtained a system of partial differential equations: 
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= d l h;-h;dM^A^lP = K*-Kl^V\^h\y 

= 3h* 5 /2h 5 -hl/h 4 
? 0, h* 5 ^ 0, 



(52) 

(53) 
(54) 
(55) 

(56) 



cases with vanishing h\ and/or h\ should be analyzed additionally. 

The system of second order nonlinear partial differential equations (|52jl - 
()55|) can be solved in general form if there are given certain values of func- 
tions g 2 (x 2 ,x 3 ) (or, inversely, g 3 (x 2 ,x 3 )), /i 4 (x\t>) (or, inversely, h 5 (x\v)), 
oj(x\v) and of sources T 2 (x 2 ,x 3 ,v) and T 4 (x 2 ,x 3 ). 

We outline the main steps of constructing exact solutions and for the 
case T 4 = when the equation (|5^|l . equivalently, (|2Tj) . is solved by the 
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h-components of d-metric g\ = 1, (72 = — 1 and 93 = ~ (% 2 ) ■ For ^2 = 0, 
the equation (J53J) relates two functions h 4 (x\v) and h 5 (x\v) following two 
possibilities: 

• a) to compute 

V\h\ = h m (x*) + h m (x 1 ) J ^\h A {x\v) \dv, hi (x\v) ^ 0; 

= h s[1] (x l ) + h m (x*) v, h\ (x\v) = 0, (57) 

for some functions ^5(1,2] stated by boundary conditions; 
b) or, inversely, to compute h 4 for a given h 5 (x\ v) ,h%^ 0, 

= h [0] (x l ) (y/\h 5 (^vfi)*, (58) 

with h]n] (x l ) given by boundary conditions. It is convenient two con- 
sider the parametrization 

M = hf 0] (x l ) [r {x\v)] 2 and \h s (x\v) | = (/ (x\v) + f (x*)) 2 

(59) 

solving (|3Sjl. We note that the sourceless equation (J53|) is satisfied 
by arbitrary pairs of coefficients h 4 (x l ,v) and h^o] {x % ) . Solutions with 
Y 2 7^ can be found by ansatz of type 

h 5 [T 2 ] = h,h 4 [T 2 ] = ft (x\v) h 4 , (60) 

where h 4 and h 5 are related by formula (|H7jl . or (JHHj) . Substituting (jHTljl . 
we obtain 

ft (a?*,v) = ft[o] (^) - y T 2 (x 2 ,x 3 ,v)-^-dv, (61) 

where ft[o] (a; 1 ) are arbitrary functions. 

The exact solutions of (JH3J) for /3 7^ are defined from an algebraic equa- 
tion, WiP+ati = 0, where the coefficients j3 and a» are computed as in formulas 
(jo^j) by using the solutions for (J52j) and (J53j) . The general solution is 

= «9 fc ln[V|M 5 |/|^|]/^ MVMsl/l^l], (62) 

with = d/dv and ft,g 7^ 0. If h* 5 = 0, or even h* 5 but /3 = 0, the coeffi- 
cients Wk could be arbitrary functions on (x l , v) . For the vacuum Einstein 
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equations this is a degenerated case imposing the the compatibility condi- 
tions (3 = oti = 0, which are satisfied, for instance, if the h± and h§ are related 
as in the formula (J58|) but with h[ ] (x*) = const. 

Having defined and h& and computed 7 from ()56j) we can solve the 
equation by integrating on variable 'V the equation n** +772* = 0. The 
exact solution is 



for some functions nk[i,2] stated by boundary conditions. 

Summarizing the results for the nondegenerated cases when hi ^ and 
hi ^ and (for simplicity, for a trivial conformal factor u), we derive an 
explicit result for 5D exact solutions for gravitational algebroids given by 
ansatz (ffijj). 
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